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The equations governing the evolution of non-minimally coupled scalar matter and the scale 
factor of a Robertson- Walker universe are derived from a minisuperspace action. As for the min- 
imally coupled case, it is shown that the entire semiclassical dynamics can be retrieved from the 
Wheeler-DeWitt equation via the Born-Oppenheimer reduction, which properly yields the (time- 
time component of the) covariantly conserved energy-momentum tensor of the scalar field as the 
source term for gravity. However, for a generic coupling, the expectation value of the operator which 
evolves the matter state in time is not equal to the source term in the semiclassical Einstein equa- 
tion for the scale factor of the universe and the difference between these two quantities is related to 
the squeezing and quantum fluctuations of the matter state. We also argue that matter quantum 
fluctuations become relevant in an intermediate regime between quantum gravity and semiclassical 
gravity and study several cases in detail. 
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I. INTRODUCTION 



Non-minimally coupled scalar fields have been extensively used in the literature to account for higher loop quantum 
corrections to the theory of scalar fields in the general relativity theory of gravitation (see, e.g., [jjj and Rcfs. therein). 
However, contrary to the minimally coupled case, quantization of such fields in the framework of quantum field 
theories on a classical curved background makes an ambiguity apparent. In fact, there appear a boundary term in the 
action for the scalar field which arises when using rescaled fields and whose presence changes the canonical matter 
Hamiltonian [^|. Thus, depending on whether one retains or drops such a term yields apparently different Schrodinger 
equations for the state of the scalar field ||. 

This ambiguity can be understood by noting that, classically, the above mentioned boundary term is generated 
because the rescaling of the scalar field by a function of the scale factor of the universe is a time dependent canonical 
transformation Therefore, in the quantum theory, it is associated with a (possibly unitary) transformation between 
different Fock spaces (choices of the ground state for matter). The ambiguity can then be eliminated || by introducing 
the invariant Fock space built upon invariant operators Q which are not affected by a boundary term added to the 
action. However, since the matter energy (the time-time component of the energy-momentum tensor) is related to 
the matter Hamiltonian, the interpretation in terms of particle content (or gravitational "weight" ) of a quantum state 
of matter remains an open question. 

The latter is not an academic issue, since, classically, it is the matter energy which acts as a source for gravity 
in the Einstein equations and determines the history of the universe. Hence, it is of primary concern to clarify by 
means of which quantity a quantized state of matter drives the metric at the semiclassical level. In order to do so, we 
appeal to the principle that time-reparameterization invariance is a fundamental property of any gravitational system 
which also generates the dynamics and be lifted to a quantum symmetry. Then, one expects to have an Hamiltonian 
constraint from which the Born-Oppenheimer (BO) reduction |^| allows to properly and unambiguously recover the 
semiclassical limit starting from the Wheeler-DeWitt (WDW) equation |&ffo| . 

The above scheme has proven remarkably successful for the minimally coupled scalar field in Robertson- Walker 
space-time |Q. In that case one starts with an action S = S(N, a, <p) for the three minisuperspace variables N (lapse 
function), a (scale factor) and <fi (one mode of a scalar field) which are functions of an arbitrary time t. Upon varying 
S one obtains three Euler-Lagrange equations of motion (for the details, see next Sections and, e.g., Ref. ||) 

6 S 

-jrjy = — H = (Hamiltonian constraint) (1-1) 
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(Klein-Gordon equation) . 



(1.2) 
(1.3) 



Since Eq. (1.1) is a constraint, N = N(t) is arbitrary and one can set, e.g., N(t) = 1 for fcj < t < tf (proper time 
gauge) provided the initial conditions a (U), a(U), 4>(U ) and (p(U) are such that H(U) = (an overdot denotes the 



derivative with respect to t). Then Eqs. (1.2) and (L3) will evolve the initial data to any time U <t<tf consistently, 
since 
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which proves that the Hamiltonian constraint is preserved by the dynamics. However , for practical purposes it is more 
convenient to revert the above inference and observe that the equation of motion ( |l.2| ) for a is identically satisfied 
provided the Hamiltonian constraint (1.1) is enforced at all times along with the Klein-Gordon equation (1.3), 



H(t) = 



^=0 



(U < t < t f ) 



<55 
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(1.5) 



A further simplification is then obtained at the semiclassical level, where one can show that both the Hamilto- 
nian constraint and the Klein-Gordon equation follow from the WDW equation obtained by quantizing the super- 
Hamiltonian H, 



HV= (H g + Hm) 9 = 0, 



(1.6) 

where Hq and Hm are respectively referred to as the gravitational and matter Hamiltonian. In fact, the BO decom- 
position for the wavefunction of the universe, 



*(«,<£) = ip(a) X {4>,a) , 



leads to the coupled equations 

H G +(H M ) =0 



i ft 



dt 



Hm Xs 



(Einstein-Hamilton- Jacobi equation) 



(Schrodinger equation) , 



(1.7) 



(1.8) 



(1.9) 



along wit h ex plicit conditions on the wavefunction s tp and \ f° r the validity of such an a pproximation. Wc then note 
that Eq. (1.8) is the semiclassical analogue of Eq. ( fj~l| ) and the operator Hm in Eq. (1.9) evolves in time states \s of 
t he s calar field in such a way that the expectation value of <f> over coherent states satisfies the Klein-Gordon equation 
(1.3). If one regards coherent states as being the quantum states which are clos est to classical, one can conclude that 



the entire semiclassical dynamics is encoded in the single WDW equation (1.6) and is properly retrieved by the BO 
reduction. 



We wish to remark here the fact that the d ecom posi tion (1.7) is not related to the gravitational scale of mass being 
bigger than any matter scale, therefore Eqs. (1.8) and (1.9) do not require an expansion in powers of the Planck mass 
[ p2| (see also ||, where such an expansion is shown to violate unitarity when incorrectly performed). Indeed, the 
relevant ratio is between the energy of each quantum of <p and the total energy of such quanta. This can be understood 
if one considers that the scale factor of the universe is a collective variable associated with the total energy of matter 
in space and that such an energy is presently much bigger than the energy of each of its microscopic constituents 
(described by the degree of freedom </>). Basically, this fact, the huge amount of matter particles, is the reason we 
live in a semiclassical universe p| . Consequently, one expects a failure in the semiclassical approximation which leads 



to Eqs. (1.8) and (l.£) for matter states x containing a small number of quanta. In fact, this expectation has been 
systematically verified in all the cases studied so far [OhIm . 
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When one applies the same scheme to the non-minimally coupled case, a new ambiguity arises because then there 
is no clear way of splitting the action into a matter part and a gravitational part and a preferred classical form for the 
energy-momentum tensor of the scalar field can indeed be singled out only by requiring covariant conservation Jl7| . 
Correspondingly, one could write many (classically equivalent) Hamiltonian con stra ints which, at the quantum level, 
become inequivalent "WDW" equations. Interestingly, the BO factorization (1.7) pinpoints a speci fic fo rm of (or 
operator ordering in) the WDW equation in order to ensure the existence of the Schrodingcr equati on (1.9). Further, 
the BO reduction then yields a source term for the geometry in the semiclassical Einstein equation (l.S) which can be 
easily related to the time-time component of the proper (divergenceless) energy-momentum tensor of the scalar field. 
However, such a source is not equal to the expectation value of the Hamiltonian operator in Eq. (1.9). It is only for 
the case of minimal coupling that the super-Hamiltonian is the sum of the time-time components of the (covariantly 
conserved) energy-momentum tensor and Einstein tensor and then the expectation value of the generator of the time 
evolution for the matter state is the semiclassical source of gravity. 

In the following, we shall argue that the difference between the source term in the semiclassical Einstein equation 
and the operator of time-evolution in the Schrodingcr equation is associated with the different ways matter and gravity 
are affected by quantum fluctuations around the mean value of the matter field and by the presence of a squeezing 
p8| term in the Schrodinger equation. This will lead us to define three different regimes of approximations, namely 
quantum gravity, semiclassical gravity and an "intermediate" regime in which quantum gravitational fluctuations are 
negligible but the trajectory of the (collective) gravitational degree of freedom senses matter fluctuations. Similar 
results were found previously in different approaches and, for example, the intermediate regime is called stochastic 
gravity by Hu and collaborators (for a review see e.g., [jl9| and Refs. therein). 

The approach employed in the present paper might look totally different with respect to the procedure of estimating 
the back-reaction within the framework of quantum field theory by studying perturbations of the Einstein equations 
(or in the Feynman path integral) around a given classical solution p(| . It has certainly the shortcoming that we start 
from an effective minisuperspace action in which the degrees of freedom of the system have been reduced by symmetry 
arguments prior to quantization and subsequent semiclassical approximation, rather than from the quantized set of 
equations derived from the full Einstein-Hilbert action. However, we point out that the more standard approach of 
perturbation theory also reduces the degrees of freedom of Einstein gravity plj , since it assumes the existence of a 
classical (saddle point) solution (background manifold and metric) from the onset, leaving as remnant gauge freedom 
only coordinate transformations and small diffeomorphisms of the background manifold ]2C[ |. Both approaches are 
thus questionable if one wishes to quantize Einstein gravity, but can be regarded as hopefully reliable whenever one 
aims at describing gravity in a semiclassical state as we wish to do here. 

In the following Section we start from the minisuperspace action for a mode of a non-minimally coupled massive 
scalar field in a Robertson- Walker space-time and pursue the standard canonical formalism in order to obtain an 
Hamiltonian constraint and the corresponding WDW equation. Then in Section III we apply the BO approach and 
obtain the Hamilton- Jacobi equation for the scale factor of the universe and the Schrodinger equation for the state 
of the scalar field. By making use of invariant operators we solve the Schrodinger equation and show that the source 
term in the semiclassical Einstein equation is the semiclassical extension of the time-time component of the covariantly 
conserved energy-momentum tensor. In Section ^ we specialize the results to the massive minimally coupled case, 
which we briefly review, and to the particular cases with £ = 1/6 and 1/4 which we analyze in more detail. Finally, 
in Section [v| we summarize and comment on our results. 

We follow the sign convention of Ref. Ill) and define k = 8 n G. 



II. MINISUPERSPACE ACTION 



In this Section we shall show that the classical dynamics of a non-minimally coupled scalar field in Robertson- 
Walker space-time is determined by the Hamiltonian constraint and the Klein-Gordon equation, thus generalizing the 
result (1.4) for the minimally coupled case as described in the Introduction. 

The (volume part of the) action for the non-minimally coupled real scalar field <I> = $(a;) in a generic four- 
dimensional space-time M. with metric g is given by |l|,p2| 



d 4 x- 



M 



K 



(2.1) 



where g = det g, R is the scalar curvature, /i = \j 1$ the inverse of the Compton wavelength of $ and £ a dimensionless 
parameter such that £ = corresponds to the minimal coupling and £ = 1/6 yields the conformal coupling [p3[. 
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It is possible to reduce the above action by assuming spatial homogeneity and isotropy so that M. admits a preferred 
foliation into spatial hypersurfaces of constant time t and the four- metric is given by the Robertson- Walker line element 

EH 



ds 2 



-N dt 2 



dr 2 



1 — k r 2 



sin 9 dip 2 



(2.2) 



with k = 0, ±1 for flat, positive and negative spatial curvature, 9 and <p the usual angular coordinates and r £ [0, r&) 
with r+i = 1, 7'o — r—i = +oo. The scalar curvature is then given by 



R -Jp u 



diV d 2 



aN a 2 



— + N 2 — 



We also expand the real scalar field in spatial Fourier modes, 



^pi 



(2.3) 



(2.4) 



where V = a 3 V is the spatial volume of the universe, and separate the real from the imaginary part, 



1 4) 



This decomposition yields an effective action for two minisuperspace variables a — a(t) and 
see Appendix |A|), 



a a 



TV fc a + a 2 — — 



V N 



(2.5) 

</>(i) (for the details 
(2.6) 



where v = V/n. The action S can be used to analyze (the real or imaginary part of) any of the modes with an 
effective frequency cu as given in Eq. ( A12 ). 

The above action contains both a second time derivative of a and a first time derivative of N in the same term. The 
former would cause problems with causality and requires a modification of the standard Euler-Lagrange equations 
of motion; the latter breaks the presumed time-reparameterization invariance of the system. However, we observe 
that upon integrating by parts the last term and neglecting the integrated part as dynamically irrelevant (for further 
explanation see Appendix pi), one finally obtains 



5' 



1 



Ndta* 



N 2 



2 ±2 
id (f> 



G 



N 2 
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<j> a (j> 
a N 2 



(2.7) 



in which there are no second time derivatives and dr = N dt is the proper time measure. This is the action we regard 
as properly describing the dynamics of the coupled variables a and <p. 



A. Lagrangean dynamics 



The Euler-Lagrange equations of motion following from the action S are given by 

1_ 6S_ = _H 

a 3 SN ~ a 3 



i ss _ 

a 2 da 



a 3 8(f) 
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(2.8) 



(2.9) 
(2.10) 
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where we have set N = 1 after the variation to give the expressions a simple form. This choice is consistent with the 
fact that the action (2/7) does not contain time derivatives of N and Eq. fl2.8|) is then the Hamiltonian constraint. Of 
course it must be preserved in time and, in fact, 



dH 



I a (v - £_6 2 ) (a 2 + 2 a a + k) + | a 2 a (j) 2 + u 2 6 2 ^j + a 3 <$> (j) + oj 2 <j)j 



+6£,a<jxj) (3d 2 + k) + 6£a 2 la <j) 2 + 6 a 6 + 6 a <j>) ~2uP ■ PaS 2 

-<j> — 
5a 86 



(2.11) 



vanishes identically by virtue of Eqs. (2.S) and ( 2.10 ), thus generalizing to arbitrary £ the results (|l.4| ) and (1.5) valid 
for £ = 0. 

The next step is to quantize the system and show that the WDW equ atio n encodes the entire semiclassical dynamics. 
In order to do so, one needs to consider the Hamiltonian form of Eq. (2.8). 



B. Hamiltonian dynamics 

The canonical momenta conjugated to N, a and (f> are given by 

Pn = 



(2.12) 



P a = -6a(v-^ 2 ) JL+Qta 2 ^ 



P lP = a 3 -^-+6^a 2 6 — 



N 



N 



The action (2.7) can then be written in canonical form as 



S 



dt 



P N N + P a a + P^j> -NH(P a , P h a, , 



where the super-Hamiltonian already given in Eq. (2.8) takes on the rather complicated canonical form 



H 



[aP a -6^P^>Y 



p2 

6ka (v-£6 2 ) + ^L + a 3 Lj 



2 1 6a 3 [v-£ (1 — 60 6 2 \ 
Several remarks are in order. First, the cases £ = 0, 



,3 , ,2 j2 



H = \ l-lL-Qkva+^+ a 3 uj 2 6 2 
2 6 v a a J 



and £ = 1/6, 
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(aPa-SP^Y 
6 v a 3 
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6kva+^- + a (a 2 uj 2 + k) 6 2 \ , 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



(2.17) 



(2.18) 



are clearly special since these values of £ simplify the form of the kinetic term in H , although it is only for £ = that 
the kinetic term is diagonal in the momenta ]24| . 

Second, the rr-component of the unique divergenceless energy- momentum tensor as computed according to Ref . p7| , 
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4> 2 + u 2 6 2 ) + 6£- 

2 \ / a 
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is related to the variation of the action with respect to the metric by 



(2.19) 
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It is however the (non-conserved) quantity in the right hand side (r.h.s.) of Eq. ( p. 20 ) which naturally appears inside 
H, as is apparent from Eq. (2.~ 



5 = («-f. 
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where G TT is the rr-component of the standard Einstein tensor. From Eq. ( |2.11 ) one concludes that 

4>H' , 



dH' 



(2.22) 



(2.23) 



therefore G TT = kT TTi or H' = 0, is an equivalent statement of time-reparameterization invariance (assuming v ^ 

£<2> 2 ). 

In order to lift such a symmetry to the quantum level and obtain the WDW equation one might choose either H 
or H 1 (or any other classically equivalent expression), thus obtaining quantum mechanically inequivalent "WDW" 
equations. However, the existence of the semiclassical limit via the BO reduction places some restrictions on the form 
of the Hamiltonian constraint. In particular, in order to recover a Schrodinger equation from the WD W equ ation, it 
is necessary that the coefficient of P 2 does not depend on the matter degree of freedom cf> (see Section III B for more 
details). This singles out the preferred (classical) expression 



H = W(.H 



6 a 3 



P 2 



-f +a 3 u> 2 (j) 2 -6fca (v-£(i> 2 ) 



(2.24) 



where = W ( {<j>) = v - £ (1 - 6f) 4> 2 . 

We proceed to analyze the quantum version of the Hamiltonian constraint with the modified super-Hamiltonian 
( 2.24 ) in the next Section where we shall also show that the BO approach guarantees that the metric be driven by 
the proper (conserved) semiclassical (rr-component of the) energy momentum tensor for the scalar field. 



III. SEMICLASSICAL EQUATIONS 



The quantization of the Hamiltonian constraint H = 0, with H as given in Eq. (2.24), is formally achieved by 
introducing the operators a, <j>, P a and P^ which yields the WDW equation 

i?(F Q ,P ,a,0)* = O , (3.1) 
where ^ = ^(a,4>) is the wavefunction of the universe. 



When dealing with Eq. (3.1) as it stands, one has to face several formal problems. First of all one would like to 
describe the system by means of Dirac variables, but, obviously, this is not the case for a and <fi and their momenta 
which do not commute with H. This problem could be solved by trading the original canonical variables for their 
initial values (the so called perennials; for a review see Ref. p5|). Then one wishes the canonical variables map into 
Hermitian operators and, further, needs to make sense of the ordering in the kinetic term. 

Since we are interested in the semiclassical limit for the variable a p6[ , at each step we shall assume the ordering 
which best fits to the computation and define a scalar product in the variable <p at fixed a as 



G 



/+oo 
d4>**(a,<t>)V(a,<l>) , (3.2) 
-oo 

which renders the operator 

P^^-iTidtj, (3.3) 

Hermitian provided the functions \& = ^(a, <ft) are summable in <p G R for any allowed (fixed) value of a. Analogously 
we define 

P a = -ihd a (3.4) 

and a and <f> as multiplicative operators. Since the range of a is R + one should carefully discuss the dependence of W 
on a, however, again the fact we want to recover the semiclassical limit will ease this issue because, strictly speaking, 
there is only one allowed value of a at a given time along a (semi) classical trajectory. Hence, it will practically be 
sufficient to consider small intervals of R + (at a time) . 

To make all the above more concrete, we consider the BO factorization (1.7) for the total wavefunction into matter 
and gravity parts and assume the matter functions are normalized in the induced scalar product 

( I V) )(a) r + °° 

Then we factor out the geometrical phase by defining 

, , (x\Pa\x)da> 
X = e n J X 

- + l T I (x\Pa\x)da' 

ijj = e h J V , (3-6) 

so that (x |P»| X) =0. 



A. Gravitational equation 



Upon substituting the above definitions into the WDW equation ( |3.l|) and contracting over ( % | on the left then 
yields an equation for the gravitational part 

o i ~T~ + 2 £ \<P p <t>) — ■ 

2 6a a z a 



= \ ^ ~H (4>P*Pa)\t> = &i 9) t>, (3.7) 

where ( O ) = (x \0\ x) f° r an y operator O and = v — £ ( 1 — 6 £) </> 2 . 

The term A { a 9) in the r.h.s. of Eq. ( |3.7[ ) can be regarded as describing fluctuations of the gravitational degree of 
freedom in the following sense. Upon assuming the space of the matter functions x admits a complete orthonormal 
basis | n ) and inserting an identity one obtains 

(X imiA) ( p a) q \ X) = E (X l/(M,-P*)l n){n \{P a f\ x) 

71 

~2J (n,a + q5a n \ x,a) , (3.8) 

n 

where / is any function of a, (f> and P^, q is a (positive) integer and Sa n ~ (x |/(a, 0, ^>)| n). Therefore, when 

is not small [with respect to the left hand side (l.h.s.) of Eq. (p?7j)] the system has a non-negligible probability of 

spreading over states n(<p;a') with a' ^ a. 
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On the other hand, when Ai 9 ' is negligible p^] one can assume x is peaked on a given trajectory a — a(r) and is 
well approximated by the WKB form 



IpWKB 



P a (a) da 



(3.9) 



where P a = Po,{t) is the momentum along the s emiclassical trajectory and the time r is correspondingly defined 
according to the semiclassical version of Eq. ( 2.13 ), 



P, 



-6 a ( W £ ) d 



For 



Pai'WKB = Pa IpWKB , 

one obtains the semiclassical Hamilton- Jacobi equation for a, 



(3.10) 



(3.11) 



P 2 
6 a 



6£ 



+ a 3 ( 2 ) = , 



(3.12) 



or, after substituting for P a from Eq. (3.1C) 



3 n 



a 3 W 2 



:w e : 



-3r 



: p2 



(3.13) 



This expression shows a deep entanglement between the two degrees of freedom of the system, so tha t it is not possible 
to distinguish a gravitational Hamiltonian from a matter Hamiltonian uniquely, as was done in Eq. (Q, when f ^ 0. 

The meaning of in the r.h.s. of Eq. ( 3.13 ) is that it describes quantum fluctuations of the matter field. In fact, 
such a term corresponds to higher 7i-order terms in the expansion of x around a classical state, for which 



P 



0(h) 



(3.14) 



thus A0 measures departures from classicality of the scalar field and could also be associated with fluctuations of the 
"effective" Newton constant, which would explain why A^ ~ £ 2 These fluctuations are formally different from 

the one s de scribed by A a 9 ^ and associated with the classicality of the gravitational degree of freedom which appear 
in Eq. (3.7). However, it is also expected that both kinds of fluctuations must be consistently small in a semiclassical 
regime. In fact, terms like A^ are usually absent in the standard expressions for the source in the semiclassical 
Einstein equations as obtained in q uantum field theory in Robertson- Walker space p]Jl9|. 
When A^ is negligible, Eq. (|3.13|) simplifies to 



3 a 



[Ws)a 2 + k (v-£($ 2 )) 



Pi 



, 3 2 

+ a oj 



(3.15) 



which then becomes equal to the classical Hamiltonian constraint (2.S ) if we replace ( (f> 2 ) with 2 and ( P 2 ) with P 2 . 
This shows that the BO approach yields the correct classical limit for gravity. 

At this point it might help in further understanding the three equations (3/7), ( 3.13 ) and ( 3.1 5| ) derived so far 
to recall the three regimes of approximation which were mentioned in the Introduction and give for them explicit 
definitions in terms of the relevant quantities: 
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1. when a! 9 ' and are not negligible [with respect to the corresponding l.h.s.s in Eqs. (3.7) and ( |3 .13 )], one 
is at the border of quantum gravity in which gravitational fluctuations start to be si gnifi cant. All the present 
treatment of the coupled matter-gravity system, and including the starting action (2/7), is likely to lose its 
meaning when A^ 9 ^ and A^ are very big, but one can still try to use Eq. (3.7) whenever such terms are not 
overwhelming; 



2. when both A^ and are negligible, one is at the opposite limit of semiclassical gravity described by Eq. (3.15). 
This provides a good picture of present day universe and already contains the description of important effects 
related to the quantum nature of matter such as the production of particles induced by the evolution of the 
scale factor a — a(t) [PJl5|l; 

3. when is negligible, but A^ is not, one is in the "intermediate" re gime where matter fluctuations play a 
significant role in driving the evolution of the metric according to Eq. ( p,13| ). One might expect that this is a 
fairly good settin g fo r the description of an early epoch in the history of the universe, e.g., near the time of the 
onset of inflation jig] . 

We also point out that, in the above scheme, one could include within matter fluctuations the effect of higher WKB 
orders in the expression for ip (representing "collective gravitational fluctuations"), which might play a significant role 
at both stages | and | f||. 

We finally note that, so far, we have not chosen a specific ordering between <f> and to keep the discussion as 
general as possible. 



B. Matter equation 



Let us now turn to the equation for the matter state. Upon subtracting Eq. (3/7) from Eq. (|3.l| ) and dividing by 
tp — ipwKB one obtains 



Pa 
1 

' 2 



3 a a A a z 

j2 I d2 



+ H' S -{H' S ) x 



P 2 - (P 2 ) £ / 

" -2^ (oP P a -{oP P a 



3 a 



where we have also used Eq. ( 3.11 ) and 

+ a 3 LU 2 4> 2 ~6ka (v-C4> 2 ) 



P 2 



A 2 p2 



3£ 2 



3fcflu 



From Eq. ( 3.10 ) one obtains 



d d 
P a P a ~ ifra— = i?i— , 
oa or 



(3.16) 



(3.17) 



(3.18) 



and this explains our previous statement that one cannot allow for a factor /(</>) multiplying P 2 in the super- 
Hamiltonian. In that case one would have 



f(<j>)P a PaX~f(<f>) 



dr 



(3.19) 



Then, the only way of getting rid of the operator f(4>) and obtain a Schrodinger equation is to assume the ordering 



d 



am 



(3.20) 



and write 



(3.21) 







However, one would eventually find that this procedure is equivale nt to the (simpler) choice of multiplying H by 
1/ f(4>) from the onset, as was indeed do ne by introducing H in Eq. ( 2.24 ). 
Upon substituting in for P a from Eq. ( |3.1C ) yields 



OT 



'He 



(3.22) 



where 




6k( 



3fc< 



a 



and 



W^a 3 



pi _ I p2 

12a 



represents the effect of gravitational fluctuations on the evolution of matter as A^ does for the gravitational state 
Finally, by neglecting A« and rescaling the matter wavefunction, 

i 

Xs = e h 

one obtains the Schrodinger equation 



(3.23) 



(3.24) 



dT(H s ) -- I dT(Hs-ihd T ) 
X = e 11 



X - 



(3.25) 



i n -— - = H s Xs 
or 



(3.26) 



One can make use of Eq. ( 3.1 3| ) to substitute for a in the operator H$ which appears in the r.h.s. above and generates 
the evolution of the matter state along the semiclassical trajectory a — a(r), thus obtaining a fairly complicated 
expression which does not naively compare to the matter source that appears in the semiclassical Einstein equation 

For a generic value of £ it looks hopeless to find invariant operators H for the complicated Hamiltonian Hs, but 
it can be done at least wh en ter ms of order v^ 1 are negligible (along with ~ A^; see case @ in the scheme of 



approximations of Section [II A). Then 



and Eq. ( 3.23| ) reduces to 
H s 



Pi 



+ a (a 2 oj 2 + 6£k) 2 -6^ ^P^ + P^^) 



(3.27) 



(3.28) 



where we have also symmetrized the product <p P^. Such a term is related to the squeezing of the matter state (see, 
e.g. the review [Q) and it is remarkable that it appears in the operator that evolves the state of the field <fi. In fact, 
upon quantizing the scalar field in a Robertson- Walker universe, one obtains a squeezing term in the Hamiltonian for 
the rescaled field ( = acf>, which is however absent in the Hamiltonian for the unsealed (f> || . Since the squeezing is 
also related to the decoherence of classical solutions |3(],|l9| , this might have interesting consequences for the onset of 
a classical universe. 

It is important to note that neglecting terms of order u _1 = k/V is tantamount to perform an expansion in the 
Newton constant and ne glect terms of order k and higher. In fact, a (infinite) factor of V has been absorbed in the 

_1 henceforth. As explained in 



definition of <fi 2 [see Eq. (2.4)], therefore one can formally set V = 1 and take v 



Ref. |8J, the expansion in k must be performed carefully in order to preserve unitarity and, indeed, in the present 
paper we expand after completing the BO red uction which yields the Eqs. (3.13) and ( p. 23] ). Had we expanded and 
truncated the original super-Hamiltonian (|2.8|) before applying the BO approach might have led us to an unphysical 
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picture in which, e.g., one neglects the matter source (of order re ) with respect to the gravitational part (of order 

For the above approximate Hamiltonian H$ one finds that the annihilation and creation operators a^ and are 
given by 



0£ 



a 3 W£ 
2h 



i 



where the effective frequency is 



wf = w 2 + 6£ 



fc St" 



The latter quantity must be strictly positive, that is |3lJ 

■2 



36£ 2 — < w 2 + 6£ — 



(3.29) 



(3.30) 



(3.31) 



which is always satisfied for a minimally coupled massive (p ^ 0) scalar field, but places restrictions on the (Hubble) 
ratio h = a/a for all other values of the parameters p, P • P and k. In particular, the case u> = k = is excluded from 
the present analysis. The Hamiltonian can then be written as 



(3.32) 



however, this does not imply that the model contains states with a fixed number of quanta. 

In fact, exact solutions of Eq. ( 3.26 ) are constructed from the invariant operators and bl pH , 

1 f 4> 



b e = 



r Th Pi 



according to 



where 



\n,r)=e 



p$ P4, - (a 3 p ? + 6 £ a 2 a p^) 



dr' s ( n, t' I Hs — ihd T \ n,r' ] 



(3.33) 



(3.34) 



(3.35) 



and we remark that the phase in Eq. ( |3.34 ) is the same that relates Xs to (the time-independent) \ in Eq. (3.25). 
The function p^ = Pi(r) is a solution of the equation 



Pi + 3 - p 5 



£ [R + 6 (1-60 



(3.36) 



with R the scalar curvature (2.3) in the gauge N = 1. We observe that setting the r.h.s. of the above equation to 
zero yields the classical Klein-Gordon equation ( [2.10 ) only for the two cases £ = 0, 1/6. 
By introducing the new variable 



_ 3/2 

<7£ = a ' p 6 



(3.37) 



one simplifies Eq. ( |3.36| ) to the form 



h + ^t^i = -3 



(3.38) 
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where 



^ 2 =^ 2 + 6 ^ + 7 I 16 ? (1-30-1] -2 + (4^-1) - 



(3.39) 



We then note that another special case is obtained for £ = 1/4 in flat space, since Sly 4 — lu for k = 0. This case is 
of interest also because every solution 77 of the (free) Dirac equation in curved space-time satisfies the Klein-Gordon 
equation 



(3.40) 



where = 1/4 [|l[|32| . Therefore, one can consider the scalar field with £ = 1/4 as the analogue of a minimally 
coupled fermion field. 



Provided is negligible, a solution to Eq. ( p.38|) is given by 



<7£ ~ SI ] 



-1/2 



(3.41) 



One then finds that s (n,T \N\ n,r) s = s (n,T \a^a^ \ ti,t) s is generally a time-dependent quantity. Also, in general 
there exist suitable initial conditions (at r = Tj) such that 



1 



3d(r ) 



k( T o) = ad T o) ) 



(3.42) 



2a(r ) 



and the state | ti,t) then describes exactly n (Hamiltonian) quanta at r = tq. In the following we shall often find it 
convenient to consider tq 3> Tj corresponding to a period of (almost) adiabatic expansion with a -C 1. 

We conclude this part by recalling that it is usually stated that the correct way of counting particles in quantum 
field theory on curved backgrounds is not by means of the number operator N ~ Hs but by introducing a (localized) 
detector coupled to the matter field [Q. In the present context, one might also dispute that ther e is t he alternative 
option of counting the number of particles by making use of their "weight" as it appears in Eq. ( |3.13 ), since this is 
what drives the observable a. We shall further investigate this point in the next Sections. 



C. Coupled dynamics 



In the same (sm all k ) approximation ( 3.27 ) which led us to the Schrodinger operator Hs in Eq. ( 3.28 ), the gravi- 
tational equation (|3.15|) becomes 



a 2 k \ (H M ] 



(3.43) 



where we have neglected ~ v 2 and (v ~ k 1 ) 



H M = - 



P 2 



-^■ + a (a 2 w 2 + 6£fc) 



(3.44) 



It then follows that the expectation value H s = s (n, r \Hg \n,r) can be expressed in terms of the function as 



TT h ( 1 

H. = ~ \n + - 
s 2 V 2 



k 
a 2 



^ + 6£-+( i -36^) - 



2 -20^ ■ 1 I 

°"? + °e - 3 - (T t <J t + ^2 J ' 



(3.45) 



and, after subtracting the squeezing term, one finally obtains the "weight" H m = s {n,T \Hm I n i T ) s as 



n ( 1 
H m = -\n+- 



k (% 
a 5 



w'+Bf^ + lj+lB* (26-1)) - 2 



<j\ + a 2 + 3 (4£ - 1) - a* &c + i I . (3.46) 
55 a a? 
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At least for a large, one expects that terms proportional to a and higher time derivatives of a become small [with 
respect to u> or k/a, see condition ( 3. 31 )]. In this approximation one then has, to next to leading order, 



•6£- 



-1/4 



H- 2 



-5/4 



3 1 
+ 16- 

' PP 
~2^~ ' 



16^(1-30 d 2 



U3' 



•6£- 



3 

8 lu'- 



(3.47) 



The next step is to substitute H m in Eq. (3.46) with the above expression for er^ into Eq. (3.43) which yields a master 
equation for the scale factor of the universe a — a(r). 

This master equation can then be integrated (at least numerically) for different choices of the parameters. The 
latter must be so chosen that the following conditions hold: 

1. {a/a) 2 and a/a must b e sm all with respect to lo 2 + k/a 2 , as is required by the condition (3.31), and the 
approximate expression ( 3.47 ) can then be employed; 



2. the number of invariant quanta n 3> 1 in order for Aa and Aa to be negligible |l3-16 ; 

3. the effect ive volume of the universe va 3 /n 3> 1/w so that A^ is negligible and Hs can be approximated by 
Eq. ( 3.28 ). This last condition can also be written as 



n < 



(3.48) 



which is like a dilute gas approximation, and £ p = Vhn is the Planck length. 

In the following Section we shall try and relax the condition || to estimate the effects of matter fluctuations. 

As we mentioned at the end of Section [II B , one could think of using the quantity H m to measure the actual energy 
of matter in the universe. Thus, one is led to identify a 3 (T T T ) = V (T T T ) = —H m and impose the conservation law 



V-(T)=0, 



(3.49) 



which is consistent with the fact that the l.h.s. of Eq. (3^43) is recognized as the (rr-component of the) stan dard 
Einstein tensor for the Robertson- Walker metric [see Eq. (|2.22| )1. As anticipated at the end of Section II B , this 
shows that the BO reduction automatically provides the correct expression for the energy-momentum tensor in the 
semiclassical Einstein equations. 

Since the scalar field is (classically) equivalent to a perfec t flu id ]ll]jL7j |, we can assume the preferred foliation of 
the space-time M. in which the four-metric takes the form ( |2.2|) corresponds to the frame comoving with this fluid 
of (gravitating) energy H m (see note p6||). The (semiclassical) energy- momentum tensor of the scalar field in the 
comoving frame can be written as pT| 



[T) =diag 



H„, 



,P,P,P 



(3.50) 



where p is the pressure. The spatial components of Eq. (3.49) then imp ly that p = p(r), in agreement with the 
hypothesis of homogeneity and isotropy, and the r-component of Eq. ( |3.49 ) yields the expected relation 



P 



6H m 
SV 



3 d a 2 



(3.51) 



which can be used to determine the pressure once H m and a have been obtained. Here we only observe that one has 
dust (p = 0) whenever H m is constant. 



IV. SPECIAL CASES 



It is clear from the previous analysis that the cases £ = 0, 1/6 are particularly simple for a variety of reasons, 
including the fact that Wq = W\i§ — v. Thus, we now review the minimally coupled case and study in detail the 
conformally coupled case. We shall also consider the case with £ = 1/4 (the analogue of a minimally coupled fermion 



field) for which the invariant structure is particularly simple in the limit (3.27) when P = k = 
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A. Minimal coupling 



For £ = the separation between matter and gravity is clear from the outset, since the semiclassical equation (3.15) 
for a is given by 



PI 



3 a (a 2 + k) = - 



where the rr-component of the energy-momentum tensor, 

1 ((Pi" 



+ a 3 u? 



a 6 (T TT \ 



(4.1) 



(4.2) 



equals the expression that is obtained by quanti zing the field <p on the Robertson- Walker background and taking for 
the energy density the expression given in Eq. ( 2.19| ) for £ = 0. This shows the equivalence of the BO approach for 
the minimally coupled scalar field to the computations performed in the more common framework of quantum field 
theory on curved backgrounds. 

O ne als o has the identity k Hs = Hm and the invariant annihilation operator, for the cases when u> ^ [see 
Eq. (3.31)], reduces to the simple form 



'2h 



Po 



+ i I po P,p - a pQi 



with the function po = Po(t) determined by the equation 

„ d . 2 1 
po + 3- Po+lu po = -t— 3 . 
a a° pQ 

This yields the exact (invariant) Fock space of states | n) with which one can also build coherent states, 

b \a) = a\ a) , 

such that the expectation value 

C (r) = s (a,T \ a,r) s , 



satisfies the classical Klein-Gordon equation ( |2 . 1 Oj ) for £ = |15j 



+ 3 - 4> c + w <j) c = 
a 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



This is the last step required to show that the semiclassical dynamics of a and <fr can be retrieved from the WDW 
equation alone when the scalar field is minimally coupled. 



For the homogeneous mode in flat space, P = k = 0, the approximation (3.47) yields 

9 el a 2 ' 



H rn He 



1 



1 



(4.8) 



where = H p is the (inertial) mass of one scalar quantum. A part from the value of the numerical factor mult iplying 
the second term inside the square brackets, H s is essentially the same as the analogous quantity computed in IV B 1 , 
we therefore do not analyze the case £ = any further and refer the reader to Ref. Jll| for its application to chaotic 
inflation. 
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B. Conformal coupling 



For £ = 1/6 one has again a considerable simplification in the semiclassical equation (3.15) for a, 

'{Pi: 



3 a (d 2 + k) - - 



+ a (a 



2 L0 2 



k) (0 2 



12 a 3 



b 2p2 



>P<h 



(4.9) 



but matter fluctuations do not disappear (A^ ^ 0). Correspondingly, the Hamiltonian which evolves the matter states 
is given by 



4.-1 



6 u a 3 



(4.10) 



and by ( 3.28| ) (with £ = 1/6) when terms of order v -1 are negl igible . In the following we shall assume such an 
approximation and use the expressions given in Sections [II B and III C for the invariant Fock space with £ = 1/6. 
We can then study the evolution of the scale factor of the universe corresponding to a matter cont ent g iven by 



a/6- 



For this we can easily estimate H s and H m by employing the expression for <J\i% given by Eq. ( 3.47 ) 



1. Homogeneous mode in flat space 



For P — k — one has 



H m ^nrrid, 1 



H s c^nm^ 1 



3g a 2 
8a 2 

5£^\ 
8a 2 



(4.11) 
(4.12) 



Upon substituting (4.11) into Eq. (4.9) we get the trajectory a = a m (r) as a solution of the master equation 



2 _ 1% 8n£ 4> a 2 
6a ~^8a 3 



n£ 2 p £ 4 



Had we used instead fl4.12|) we would have got a different trajectory a = a s (r) which solves 



3d 2 



4n<^, a 2 



£4, 4a 3 -3n^^ 



(4.13) 



(4.14) 



There is then an obvious difference between a m and a s , that is the velocity of the former is always finite for a > 0, 
while the velocity of the latter diverges for a hnite positive value of a [after the dilute gas approximation ( 3.4S ) has 
broken down]. 

In order to show the difference explicitly, we give a hrst example in Figs, [l] and ||, where we consider a couple 
of solutions of the above Eqs. ( ]4.13 ) and (4.14) with 1$ = nl^i^ — 1 (in natural units h = k = 1). Further, we 
have chosen a(0) = 1 for both equations so as to avoid the singularity in d s and compare two trajectories starting 
at the same value. In Fig. [I] we show the trajectories a m and a s along with the corresponding Hubble coefficients 
h m — a m /a m and h s = a s /a s and in Fig. || we plot H m [a m ], H s [a m ] and H* = H s [a s ]. It is interesting to note that, 
although the number of invariant quanta, n, remains constant in time for the exact solution | n) 1 ^ e , the number of 

quanta N ~ H s as computed from the expectation value of the Hamiltonian H$ decreases and, on the contrary, the 



"weight" H m of the state 



'1/6 



increases. The quantity represents what the "weight" of the state would be 



were the squeezing factor totally absent. At late times (r — > 00) terms proportional to h m and its derivatives vanish 
(adiabatic limit) and the three quantities converge to the same valu e. Th is is well suited if one aims to study the 
evolution of the universe assuming to know its present state [see Eq. fl3.42 )]. 

In Eq. ( 4.1 3| ) there are two parameters one can vary, that is £^ and n. In Fig. || we show the effect of taking 
X/lcf, = /i = 1,10,100 (with n£ 2 £<p = 1) on both a m and a s with a(0) = 1 and in Fig. || the effect of changing 
n = 1, 10, 100 (with fi = 1) on a m and a s with a(0) = 5. In Fig. |] and ^ we plot the corresponding Hubble coefficients 
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h m . In particular one can see from Figs. || and || that both the scale factor and the Hubble coefficient scale with a 
positive power of n, as one expected from the fact that H n 



Since Eq. (4.13) does not forbid a m (r) to approach zero, we also plot in Fig. |^ a trajectory with /i = nl 2 lj, = 1 



which starts at a m (0) = 0.2 and the corresponding Hubble coefficient h m . In the same graph we also show that 
t erms proportional to a 2 /a 2 and (a/a) (a 2 /a 2 ) are negligible for the trajectory a m as is required by the approximation 
( |3.47 ). Fig. H reproduces the behaviour of H m and H s . In particular, it is now apparent that H m st arts out at about 
zero and does not diverge for a m (0) — ► 0. However, in this case nlpl$ ~ a 3 (r) and the condition ( 3.4§| ) is violated. 
Thus, we have also plotted 



i\n 2 
24 a 3 



(4.15) 



to compare its relevance with respect to H m . 

Since dominates at small a, we have computed a corrected trajectory a c = a c (r) with a(0) = 0.4 which we plot 
in Fig. |^ together with the corresponding Hubble coefficient h c and a trajectory a m with the same initial condition. 
Finally, in Fig. we display the behaviour of 



H r = H„ 



(4.16) 



together with H m and H s for the trajectory a c 



2. Massless modes 



This is also a remarkable case, corre spond ing to what is usually considered true conformal coupling (because of 
fi = 0). Indeed, in the approximation ( 3.47 ) we find that there is no difference between the gravitational "weight" 
and the expectation value of the Hamiltonian H s , 



n / - - 

H s ~n - [P-P + k 
a \ 



(4.17) 



where P-P + k must be stri ctly positive [see Eqs. (3.31 and (3.47)]. 

Since the term in Eq. ( 4.15 ) again dominates for small a, the scale factor is determined by the master equation 



, nil ^ ^ nB \ 

a 2 = — | [P P + k v - - k 

3a 2 \ 24a 2 / 



(4.18) 



In Fig. 11 we plot the solution a c — a c (r) of Eq. ( 4.18 ), its Hubble coefficient h c and gravitational "weight" H c for 
the mode P ■ P = n = 1 in flat space, k = 0, and o c (0) = 0.21. 



C. Fermionic coupling 



For 5=1/4 and P = k = one finds that 



01/4 



1 

VJ 1 



(4.19) 



is an exact solution of Eq. (3.38) which holds in the approximation (3.27). It then follows that 

H m =H s =n m(t> , (4.20) 
which is a constant. The corresponding master equation for the scale factor, once one includes A^, is given by 



o • 2 v 
6a = n — 



Is 



24 a 3 



(4.21) 



We plot the solution a c for a(0) = 0.4, in Fig. [IJ, together with the Hubble coefficient h c and the gravitational "weight" 
H Cl and observe that the qualitative behaviour of these three quantitie s is sim ilar to the one of the analogous quantities 
for the massless conformally coupled scalar field described in Section IV B 2. 
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V. CONCLUSIONS 



In this paper we have analyzed the dynamics of a mode of a real scalar field non- minimally coupled to the Robertson- 
Walker metric. Starting from the Wheeler-DeWitt equation, we have employed the Born-Oppenheimer approach which 
has then led us to a semiclassical picture in which the state of the scalar field is evolved by a Schrodinger equation 
and the scale factor of the universe by a semiclassical Hamilton- Jacobi equation. The main result is that, for generic 
coupling, the expression for the gravitational "weight" of a matter state is not naively related to the Hamiltonian 
operator appearing in the Schrodinger equation and evolves in time differently with respect to the expectation value 
of the latter. Correspondingly, the scale factor of the universe evolves accordingly to a non-trivial master equation. 

By choosing the parameters of the model so as to obtain a Schrodinger equation for which the exact (invariant) 
Fock space can be constructed using known methods, we have studied such master equation for the cases which are 
mostly treated in the literature, that is the massive minimally coupled (£ = 0) scalar field and both massive and 
massless scalar fields with £ = 1/6. Further, we have considered the homogeneous mode of a massive scalar field with 
£ = 1/4 in flat space whose Klein-Gordon equation is formally the same as the one satisfied by minimally coupled 
Dirac fields and for which the dynamics of the scale factor shows remarkable qualitative similarities with the case of 
the massless conformally coupled (£ = 1/6) scalar field. 

For £ = 1/6 we have explicitly shown that the gravitational "weight" H m of a given massive matter state increases 
in time, at least during the early stages of the expansion. In the spirit of the principle of equivalence, according to 
which the gravitational mass of a particle equals its inertial mass, this can be considered as the signature of real 
particle production. In fact, although no local detector has been introduced, one can regard the scale factor itself as 
an observable quantity, e.g., by means of measuring the recession of galaxies, and relate the counting of particles to 
its evolution. 

We also found that H s , the expectation value of the matter Hamiltonian in the Schrodinger equation, generally 
decreases (or stays constant). Because of the different behaviours of H m and H s , one might conclude that there is 
a failure of the principle of equivalence, since one would expect that the energy by which a matter state is evolved 
in time is the same that gravitates. Were this observation proved correct, the massless conformally coupled scalar 
field and the homogeneous mode of a massive scalar field with £ = 1/4 in flat space would stand up as very peculiar, 
since for them (as for the minimally coupled scalar field) the two quantities are equal and the equivalence between 
inertial mass and gravitational mass would therefore be preserved for the massive case with £ = 1/4 (thus suggesting 
an analogous result for fermions). 

However, we point out that, while H m was shown to be the semiclassical (time-time component of the) unique 
covariantly conserved energy-momentum tensor and has naturally a physical meaning as the energy of the perfect 
fluid modeled by the scalar field, H s cannot be related to any directly measurable quantities in our treatment. Hence 
it is not clear whether H s carries any physically accessible information, although the corresponding operator Hs plays 
a fundamental role for the dynamics. In order to enlarge the number of observable quantities and make testable 
predictions one might then consider inhomogeneous fluctuations of matter fields perturbatively on the background 
determined by the master equations obtained in this paper and estimate, e.g., the effect induced on the spectrum of 
the cosmic microwave background radiation. 

We wish to conclude by mentioning that further possible extensions of the present work include a deeper analysis 
of purely quantum effec ts, su ch as those induced by t he s uperp ositi on of several matter states ]33| or the geometrical 
phase appearing in Eq. ( 3.34 ) and the r.h.s.s of Eqs. (3/7) and ( |3.16| ), and different couplings between gravity and the 



scalar field, such as those in scalar-tensor theories of gravity (for a recent review see Ref. |34|). All such extensions 
would affect the evolution of the background and, eventually, of inhomogeneous fluctuations of the matter fields. 
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APPENDIX A: ACTION FOR §p 



For £ = the scalar field potential in the action (2.1) is quadratic in $ and different modes in the sum ( [2.4| ) decouple, 
while, for £ 7^ 0, one expects the term $ 2 R induces interactions via "graviton exchange" . In any case, gravity would 
respond to the sum of all modes, thus, for the sake of simplicity, we shall consider a scalar field containing only one 
mode of fixed wave vector p , 
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$ = —/= [cos [p ■ x) 4>\ + sin (p ■ x) 4> 2 ] , (Al) 

where <pi = <ph and 2 = <Pg- The scalar product a ■ b = 7^ a 1 with 7^ = gij/a 2 — r, 9, tp) and g as given in 
Eq. (2.2), is time- independent. One then finds 



$ 2 = — [cos 2 (p ■ x) 4>\ + sin 2 (p ■ x) 02 + sin (2p ■ x) <f>x 4>2 



(A2) 



and an analogous expression for (<9$) 2 . The integration over the spatial volume yields the following constant coeffi- 
cients (7 = det 7) 



(A3) 



V c = J d 3 x -^7 cos 2 (p ■ x) 



Since the three-metric 7 is isotropic, the direction of p cannot affect the value of the above integrals so that V c and 
V s depend at most on the modulus p ■ p. Further, homogeneity of 7 implies that 



J d 3 x 1/7 sin (2p 



■x)=0 



After recalling that 

V = J d 3 x^, 

and, setting cf>2 = 0, one then obtains an action for the real part of 



Si 



Ndta A 



V N 2 



K V 



R 



with 



2 _V S p-p V c 



J1 ~ V a 2 

and, setting <pi = 0, an action for the imaginary part, 



V 



Sn = - 



1 ft 



Ndti 



Vs 4>l 2:2 



with 



2 _ V c p-p t V s 



v» 2 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



1 ~ V a 2 

For the homogeneous mode, p — 0, one has V c = V and V s = 0, so that S2 vanishes and Si coincides with the 
expression in Eq. ( |2.6| ) with <j) = 4>i an( i w = M- 

For 0, both V c and V s are strictly positive and one can rescale the fields according to 

2 1 , (A10) 



and correspondingly define an "effective" wave vector 

p =\Tv/ ( or P ^\fv/) ' 

so that the action Si (S2) for the real (imaginary) part <fi is again equal to the expression in Eq. ( |2.6| ) with 

, P-P 



(All) 



(A12) 



This shows that the action (2.6) can be used to describe the dynamics of (the real or imaginary part of) each mode 
of the real scalar field. 
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APPENDIX B: BOUNDARY TERMS 



The procedure which leads to the action (2.7) from the one in Eq. (2.6) is the analogue (for generic £) of what is 
done in general relativity when one defines the Einstein-Hilbert action as pj 



SeH -—k 



d A x- 



M 



-gR-~ 

K 



d x K 



dM 



(Bl) 



where K is the extrinsic curvature of the border dM of the space-time manifold M. In the above, the surface integral 
evaluated on dM precisely cancels all the troublesome terms inside the volume contribution (including first time 
derivatives of the lapse function and second time derivatives of the three-metric). 

For the Robertson- Walker metric (2.2) one has that K vanishes at the time-like border r — and the only 
contribution to the surface integral comes from the hypersurfaces t = U and t = tf, 
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It therefore appears natural, for £ ^ 0, to generalize the standard prescription to 
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in order to eliminate unwanted terms from the action (|2.6| ) and obtain the form (2.7). 

Of course one could also consider other ways of proceeding and, e.g., allow for terms containing d. However, the 
requirement that time-reparamctcrization remains an invariance of the system (Pp? = 0) seems to favor the above 
procedure. 
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FIG. 7. Evolution of the scale factor a m for a m (0) = 0.2, the corresponding Hubble coefficient h m and the ratios between 
velocity and acceleration of a m and jj,. 




FIG. 8. Plot of H m , H s and |A^| = — for the trajectory in Fig. ^. 
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FIG. 9. Plot of a c and a m with a(0) = 0.4 and the corresponding Hubble coefficient h c 
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